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NUMERICAL  RESULTS  OF  THE  THEORY  OF  DIFFRACTION 
OF  RADIO  WAVES  AROUND  THE  SURFACE  OF  THE  EARTH 


P-  A.  Azrilyant,  M.  G.  Belkina 


The  monograph  contains  a  collection,  in  the  form  of 
tables  and  graphs,  of  the  results  of  the  calculation  of 
the  attenuation  factor  of  the  electromagnetic  field  of 
radio  waves  with  horizontal  and  vertical  polarizations 
in  the  case  of  their  normal  propagation  around  the  Earth. 

The  main  calculations  of  the  attenuation  factor  are 
made  in  sequence  and  encompass  the  areas  of  the  shadow  and 
semishadow . 

The  joining  of  the  diffraction  curves  with  curves 
based  on  reflection  formulas  is  investigated. 

Additional  formulas,  tables  and  graphs  are  given  which 
make  it  possible  with  the  help  of  simple  operations  to  make 
a  complete  calculation  of  the  normal  propagation  of  radio 
waves  (beginning  from  the  illuminated  area  and  ending  with 
the  area  of  deep  shadow).  An  example  is  given  of  the  appli¬ 
cation  of  the  numerical  results  of  the  work  to  the  calcula¬ 
tion  of  electric  field  intensity  curves  for  certain  wave¬ 
lengths. 

The  results  obtained  are  applicable  in  the  case  of 
horizontal  polarization  for  practically  any  waves  and  any 
electric  properties  of  the  terrestrial  surface,  and  in  the 
case  of  vertical  they  are  calculated  for  the  extreme  cases 
of  very  short  and  very  long  waves,  when  the  terrestrial 
surface  can  be  considered  approximately  as  ideally  reflecting. 

The  book  is  intended  for  radio  engineers  engaged  in  the 
calculation  of  radio-wave  propagation  and  of  radar  devices. 

1 . 


1.  Basic  Formulas 


The  present  work  deals  with  the  calculation  of  the  attenuation 
factor,  characterizing  the  diffraction  of  radio  waves  around  the 
terrestrial  surface. 

Assume  r  -  spherical  coordinates,  selected  so  that  the 

emitter  is  positioned  with  t?=0.  Then  in  an  absolute  system  of  units 

the  electromagnetic  fields  of  elementary  dipoles,  emitting  radio 

waves  above  the  surface  of  the  uniform  spherical  Earth,  according 

to  the  theory  of  diffraction  have  the  following  form  (dependence 

lot  LcU-t  V 

on  time  is  understood  in  the  form  e“  se"  *  ). 

For  a  vertical  electric  dipole 


1 

x=x  — 

c 

£»  =  - 
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For  a  vertical  magnetic  dipole  (horizontal  frame) 
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(1.2) 

Here  I  -  current  in  the  elementary  dipole  or  loop,  l  -  length 
of  dipole,  S  -  area  of  loop.  Consequently  the  electric  moment  of 
the  dipole  is  equal  to  P~TL^  >  an<*  the  nagnetic  moment  of  the 
loop  fi.  denotes  the  geometric  radius  of  the  Earth. 

For  the  horizontal  electric  or  magnetic  dipole  the  formulas 
for  the  fields  have  a  somewhat  more  complex  form. 

The  horizontal  electric  dipole 

2. 
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The  horizontal  magnetic  dipole  (vertical  frame) 
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(1.4) 

In  these  formulas  it  is  assumed  that  the  moment  of  the  horizon¬ 
tal  dipole  is  oriented  in  the  direction  C]P=0c 

The  function  V  entering  into  these  formulas  depends  on  three 
dimensionless  coordinates  X,  Vi  >Va  and  the  dimensionless  complex 
parameter  q  or  q’.  These  variables  are  connected  with  the  heights 
of  the  points  of  radiation  and  reception  h-\  and  h2,  the  distance 
between  them  on  the  Earth's  surface  s,  wavelength  A  »  conductivity 
<T  »  and  dielectric  constant  of  the  soil  6  by  the  following  correla¬ 
tions  : 
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(1.6) 


3. 
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(1.7) 


where 


+  i f|/  .-l  +  i±L 


is  the  large  parameter  of  the  problem. 

Through  CL*  the  effective  radius  of  the  Earth  is  designated, 
taking  into  account  the  change  in  the  refractive  index  of  the  terres¬ 
trial  atmosphere  in  that  case  when  it  depends  on  height  linearly 
{see  [3]).  In  practice  the  values  of  CL*  can  be  most  diverse,  from 
CL*-  CL  to  d*~  00  • 

In  the  case  of  a  uniform  atmosphere  the  effective  radius  of 
the  Earth  CL*  is  equal  to  its  geometric  radius  CL  . 

The  variables  50  (A)  and  j^(A)  are  the  standard  distance  and 
height,  depending  on  wavelength  A.  and  on  effective  radius  CL*  >  i-e-, 
on  atmospheric  conditions. 

The  function  V  (  »  9 )  is  called  the  attenuation  factor. 

For  i-t  V.  A.  Fok  [13  gave  the  following  integral  representation 

- - 


■V?J 


(1.8) 


Here  T  -  contour  in  the  plane  of  the  complex  variable  from 
e  * to  0  and  from  0  to  oO  ,  function  FCt,^  >  >  d)  for  Jfe  > 

>v,  is  determined  by  the  formula 


where 


Hi, vit-y  >•»*</, 
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(1.9) 


(1.10) 


4. 


(1.11) 


and  <o  (t)  is  the  complex  function  of  Airy,  determined  by  V.  A.  Fok 
[2]  using/  the  formula 


w(t) 


1  p 

■75  J*  * 


(1.12) 


r  -  the  same  contour  as  in  formula  (1.8). 

This  function  is  the  solution  of  the  differential  equation 


*"(<)  =  tw(t). 


(1.13) 


and  u(t)  and  “t>(t)  are  those  solutions  of  this  equation,  which  in 
the  case  of  real  values  of  t  are  real  and  imaginary  parts  cu(t)  and, 
finally , 


(1.14) 


Let  us  note  that  replacement  of  the  geometric  coordinates  and 
physical  parameters  by  those  mentioned  reduces  the  number  of  variables 
on  which  the  attenuation  factor  depends.  Here  the  obviousness  of 
the  representation  of  the  function  is  in  no  way  disrupted,  since 
the  distances  and  heights  are  measured  as  if  in  a  new  scale,  depend¬ 
ing  on  wavelength  X.  and  atmospheric  conditions  (taking  into  account 
the  effective  radius  a*) .  Transition  to  coordinates  reduced  from 
geometric  is  reduced  to  the  division  of  the  latter  by  the  horizon¬ 
tal  scale  $$(2)  and  the  vertical  scale  •  ihe  scales  Sfl(A.)  and 

(X)  for  different  values  of  effective  radius  a*  are  given  in 
Figure  1 . 

According  to  the  assumption  made  in  the  conclusion,  formulas 

(1.1) — (1.4)  are  suitable  in  the  area  of  the  shadow  (without  any  limi¬ 
tations  for  values  of  the  angle  i>-  )  and  in  the  area  of  the  semishadow, 
and  in  the  illuminated  area  they  give  correct  results  only  for  a 
sliding  incidence  of  a  ray. 

For  nonelementary  radiators  the  fields  in  the  indicated  areas 
are  expressed  with  the  help  of  formulas  which  are  similar  to  formulas 

(1.1) — (1.4).  It  is  possible  to  establish  the  form  of  these  formulas 
on  the  basis  of  the  principle  that  two  sources,  giving  in  free  space 

5. 


the  same  field  of  radiation  for  the  given  azimuth  f  in  directions, 
making  up  small  angles  with  the  horizon^,  create  in  this  azimuth  <P 
in  all  the  listed  areas  the  same  diffraction  fields.  Differences 
can  be  manifested  only  in  that  part  of  the  illuminated  area  where 
rays  fall  which  are  reflected  from  the  Earth  at  angles  which  differ 
noticeably  from  ,  but  in  this  case  formulas  (1.1)— (1.4)  in  gen¬ 
eral  are  not  applicable,  and  it  is  necessary  to  use  formulas  of  the 
theory  of  propagation  of  radio  waves  over  a  flat  Earth. 

As  an  example  of  the  action  of  the  indicated  principle  one  can 
cite  the  vertical  magnetic  and  horizontal  electric  dipoles.  If  the 
dipole  moments  are  taken  such  that 


then  in  the  azimuthal  direction  ^  near  the  horizontal  plane 

they  will  give,  as  it  is  easy  to  comprehend , the  same  field  of  radia¬ 
tion  in  free  space.  Therefore  their  diffraction  fields  on  the  hori¬ 
zon  and  beyond  it  ,  as  this  is  easy  to  check  from  formulas  (1.2) 
and  (1.3),  will  coincide. 

In  the  derivation  of  formulas  (1.1)— (1.4)  it  was  also  assumed 
that  the  refractive  index  of  the  soil  ,  j  / 

is  a  large  complex  number.  * 

Now  let  us  consider  the  case  of  incidence  of  a  flat  wave  on 
a  sphere  with  a  large  complex  refractive  index. 

Relative  to  problems  of  the  propagation  of  radio  waves  this 
case  can  be  interpreted  as  a  problem  of  diffraction  around  the  surface 
of  the  Earth  of  radio  waves  which  are  arriving  from  interplanetary 
space  (nonuniformity  of  the  atmosphere  in  the  problem  of  a  plane 
wave  we  disregard). 

Assume  that  the  field  of  the  incident  wave  is  given  by  the 
formula 


(1.15) 


(  r ,  ,  <p  -  spherical  coordinates,  and  the  direction  of  propaga¬ 
tion  of  the  incident  wave  is  Then  the  diffraction  field 

near  the  surface  of  a  sphere,  at  small  distances  from  it  in  comparison 
with  the  radius  of  the  sphere,  is  given  by  the  formulas 


6. 


!»•  (*- 
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E.-Hr=*  e  - —  IMz.y.^sine 

/TInO 

i  ■'- t‘~f)  w, (».,,,)  .,.t  (1.16) 

M  sin  *> 

IMm  (t—  -  ) 
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/  m  /Tsrr  ^ 

The  function  V-j  (Z,y,  q),  appearing  in  these  formulas,  which  we 
will  call  the  attenuation  factor  for  the  field  of  an  incident  plane 
wave,  was  introduced  by  V.  A.  Fok  and  defined  by  him  in  the  form 
of  the  integral  (see  [1]) 


tu'Q(t,y,q)dt. 


(1.17) 


Variables  Z  and  y  in  this  problem  are  determined  by  the  formulas 


k(r—a) 


ka  .7 


(fr. 


(1.18) 


(1.19) 


Thus  X  is  proportional  to  the  angular  distance  or  the  distance  on 
the  surface  of  the  Earth  from  the  point  of  observation  to  the  boun¬ 
dary  of  the  geometric  shadow  1?  =  2T  on  the  sphere,  and  in  the  area 
of  the  shadow  Z>0,  in  the  illuminated  area  r  <4.  0,  and  y  is  the 
reduced  height  of  the  point  of  observation  above  the  surface  of  the 
sphere.  Let  us  note  that  since  the  nonuniformity  of  the  atmosphere 
in  this  problem  is  not  taken  into  account,  then  (1.18)  and  (1.19) 
and  the  expressions  for  q  and  q'  (1.6),  (1.7)  include  the  geometric, 
and  not  the  effective  radius. 

Formulas  (1.16),  just  as  formulas  (1.1)— (1.4),  are  applicable 
in  the  area  of  the  shadow,  semishadow,  and  in  that  part  of  the  illumi¬ 
nated  area  where  the  glancing  angle  of  the  incident  wave  is  small. 

All  of  these  formulas  are  given  in  absolute  units. 

7. 


In  order  to  switch  from  absolute  units  to  practical  it  follows 
to  express  current  in  amperes,  length  in  meters,  and  then  the  electric 
field  is  obtained  in  V/m,  if  30  ohms  is  substituted  for  1/s  [1/c]. 

In  order  to  obtain  the  magnetic  field  in  A/m  it  is  necessary  to 
replace  1/c  by  l/HTf  . 

In  the  present  work  results  are  given  from  calculations  of  the 
attenuation  factors  q )  and  q),  entering  into 

the  formulas  for  the  components  of  the  fields  of  the  vertical  and 
horizontal  (electric  and  magnetic)  dipoles  arid  the  field  of  a  plane 
incident  wave.  As  concerns  the  derivatives  of  these  functions, 
entering  into  the  same  formulas,  then  they  do  not  have  great  impor¬ 
tance,  since  the  field  components  expressed  through  them  in  the  case 
of  finite  values  of 

the  variables  Xor z  in  order  of  magnitude  are  M  times  (M  « 
large  parameter)  less  than  those  components  which  are 
expressed  through  these  functions  themselves.  Therefore  the  derivative 
attenuation  factors  q)  and  V-\(Z,y,  q)  we  did  not  calcu¬ 

late  with  the  exception  of  these  which  were  necessary  for  determina¬ 
tion  of  the  functions  themselves  for  small  values  of  heights  adduced. 

The  subsequent  paragraphs  are  devoted  to  the  investigation  of 
different  formulas  and  conditions  of  their  applicability.  Those 
interested  only  in  numerical  results  can  switch  directly  to  §6,  where 
a  description  is  given  of  the  tables  and  graphs,  and  schemes  are 
added  which  indicate  which  tables,  graphs  and  formulas  should  be 
used  in  different  areas. 

2.  Series  of  Deductions.  Dependence  on  Soil  Properties 

The  integral  representations  (1.8)  and  (1.17)  of  the  attenuation 
factors  V  and  Vi,  possessing  a  sufficient  degree  of  generality,  usually 
are  too  complex  for  direct  calculations. 

However,  these  integrals  can  be  presented  in  the  form  of  series 
of  deductions  in  points  t3 ,  being  the  roots  of  the  equation 


•'(0  —  qw(0""0. 


8. 


The  roots  ts>  corresponding  to  the  given  value  of  q,  form  a 
sequence,  increasing  in  absolute  value  and  situated  in  the  first 
quadrant  near  the  ray  arc  t  =  -~  ,  and  for  q=0  and  q=o©  -  on  the 
ray  itself. 

The  series  of  deductions  have  the  form 


=  e '  Y-gf  IzSL^zM 

®(g.  .»(o 


V,{z,y ,q)  -  i2  V^Y - — _ -  . 

(t.)  «•(*,) 


(2.2) 


(2.3) 


The  series  (2.2)  is  convergent  at  any  possible  values  of 
V,  vjz  and  can  be  considered  on  a  level  with  integral  (1.8)  by 
determination  of  the  function  q),  although  the  area  of 

practical  convergence  of  this  series  is  limited;  with  low  values 
of  X  and  large  values  of  y  it,  without  taking  the  remainder  term 
into  account,  becomes  unsuitable  for  calculations  (see  §4). 

Series  (2.3)  converges  only  when  0  and  makes  no  sense  when 
7<0,  therefore  the  determination  of  function  V i  ( JZ .,y  ,  q)  is  pos¬ 
sible  only  with  the  help  of  the  contour  integral  (1.17).  However, 
in  the  area  of  shadow  and  semishadow  these  series  converge  well. 

Let  us  consider  the  extreme  cases  of  zero  and  infinite  values 
of  q.  The  roots  corresponding  to  q=oo  we  designate  ts°,  so  that 

«(*/)- 0.  (2.4) 


and  the  roots  for  q*0  through  ts* ,  so  that 


•'(O-o. 


(2.5) 


Then,  using  equation  (2.1)  and  converting  in  series  (2.2)  and  (2.3) 
to  the  limit  qroo  ,  we  obtain 


e  4  2  /7x  **/  f.CyJ  g,  (y3) 


(2.6) 


V1(z,y,oo)-i2/r^w((<~^ 
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where  it  is  assumed  that 


e,(y)=~ 


W(t/) 


(2.8) 


For  qsO  from  (2.2)  and  (2.3)  the  immediate  result  is 

V  (x,ytyaO)  = 


1 4  **£  77-^^ 


(2.9) 


=  i  2/^V  <  -  w\fjrf,{y)' 
»1  4 


(2.10) 


where 


1 1  y)-  *(*;-?) 

My>  -5?r 


(2.11) 


Now  we  will  explain  under  what  physical  conditions  the  extreme 
cases  q=0  and  q=oo  are  realized. 

The  complex  parameters  q  and  q',  on  which  the  attenuation  fac¬ 
tors  V  and  Vi  depend,  are  connected  with  wavelength  A  and  the  soil 
parameters  V  and  t  - 

According  to  (1.7)  the  parameter  q*  is  the  product  of  two  large 
co-factors.  Therefore  it  is  always  possible  to  consider 


(2.12) 

As  concerns  the  parameter  q,  then  it  can  take  the  most  diverse 
complex  values,  which  is  clear,  for  example,  from  Figure  1,  where 
the  dependence  of  absolute  magnitude  of  q  on  wavelength  is  shown 


Figure  1.  Dependence  of  Iql  on  wavelength  ^  for  sea  water  and 
wet  soil. 

Key:  (1)  Sea  water;  (2)  Wet  soil. 


for  sea  water  and  wet  soil. 

In  [6]  an  investigation  was  made  of  the  dependence  of  attenua¬ 
tion  factors  V  and  Vi  on  the  complex-valued  parameter  on  the  surface 
of  the  Earth.  In  this  work  the  complex-valued  parameter  q  was  ex¬ 
pressed  through  the  material  parameters  ot  and  lgn  in  the  following 
manner 

’in* 


(2.13) 

The  calculations  were  made  for  values  of  o<=0,  0.01,  0.02,  0.03 
and  values  of  lgn,  changing  from  -  1  to  3»  i.e.,  for  absolute  values 
of  q,  lying  approximately  within  the  limits  from  Iql/*' 0.1  to 
Iql  /vlOO  (Figure  1).  Here  it  turned  out  that  the  dependence  of 
V  and  Vi  on  the  parameter  q  is  very  significant  when  Iql  1 ,  where 
the  moduli  of  the  functions  in  the  case  of  a  fixed  %  or  X  and  with 
an  increase  of  lgn  decrease  sharply,  and  the  phases  have  a  maximum 
(see  [6],  drawings  4-19). 
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Figure  2.  Dependence  of  attenuation  factors  on  the  complex-valued 
parameter  q  in  the  case  of  fixed  values  of  X  or  Z  and  null  values 
of  heights. 


However,  for  |q|2^1  it  is  possible  to  use  the  values  of  functions 
V  and  V-]  for  q=  oo ,  and  when  |ql^<1  their  values  for  q=0. 

Figure  2  illustrates  this  circumstance  for  heights  of  source 
and  point  of  observation  equal  to  zero.  Plotted  on  it  in  a  loga¬ 
rithmic  scale  are  the  moduli  of  functions  V  and  Vi  in  the  case  of 
fixed  values  of  X  or  2  and  the  values  of  these  functions  for  q=0. 

(On  the  axis  of  abscissas  here  the  values  of  lgn  and  their  corres¬ 
ponding  values  of  Iql  are  plotted.) 

In  Figure  3  the  moduli  of  the  derivatives 


d*V(x,0,0,g) 

dyioy. 


<7*V(jc,0,0,  q) 
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Figure  3.  Derivatives  of  attenuation  factors  in  dependence  on  the 
complex-valued  parameter  q  with  fixed  values  of  ?r  or  ;?  and  zero 
values  of  heights. 


dV^z.Q,  q) 

dy 


=  —  qVi(zfi,q) 


are  compared  for  large  values  of  |q|  (lgn)  with  their  values  when 
q=e©.  These  derivatives  on  the  strength  of  the  conversion  of  the 
functions  themselves  when  q-^eo  to  zero  determine  them  in  the  case 
of  small  values  of  y  and  large  Iqj. 

Thus  we  see  that  the  consideration  of  the  extreme  values  of 
q=0  and  q=oo  also  makes  it  possible  to  encompass  cases  when 
lq|<0.4  and  I  q  I  >40.  Furthermore,  on  the  strength  of  (2.12)  the 
case  q=  oo  makes  it  possible  to  obtain  the  main  components  of  the 
fields  of  the  horizontal  electric  dipole. 

Let  us  note  also  that  the  case  of  an  ideally  conducting  spherical 
surface  corresponds  namely  to  the  extreme  values  of  q=0  (for  vertical 
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polarization)  and  q=00(for  horizontal). 

For  the  calculation  of  attenuation  factors  in  the  case  of  arbi¬ 
trary  values  of  the  complex-valued  parameter  q  in  the  area  of  the 
semishadow,  where  at  least  five  members  of  the  series  (2.2)  and  (2.3) 
have  to  be  taken  into  account,  a  colossal  volume  of  operations  is 
required,  which  we  were  not  able  to  carry  out. 

On  the  strength  of  the  calculation  difficulties  indicated  the 
data  available  in  the  literature  are  still  very  far  from  an  exhaus¬ 
tive  numerical  investigation  of  the  diffraction  field  for  any  soils. 

Thus  in  [7]  results  are  given  (in  the  form  of  graphs)  for 
different  q  and  heights  which  are  not  equ'al  to  zero,  but  for  the 
area  of  deep  shadow,  where  it  is  possible  to  be  limited  to  the  first 
term  of  the  series  of  deductions  and  only  for  attenuation  factor 
V.  Similar  calculations  were  made  much  earlier  by  B.  A.  Vvedenskiy 
and  their  results  were  drawn  up  in  the  form  of  a  nomogram  [8]. 

In  the  book  by  Bremmer  [9]  results  are  given  in  the  area  of 
seraishadow  for  two  specific  soils  (two  values  of  £  and  C),  and  either 
for  heights  equal  to  zero  (on  various  X  )  or  for  several  heights 
not  equal  to  zero,  but  only  for  two  fixed  values  of  X  .  Here  also 
only  function  V  is  given. 

In  work  [6]  results  are  given  from  calculations  on  two  members 
of  a  series  of  deductions  for  heights  equal  to  zero,  but  for  a  whole 
number  of  soils  and  for  wavelengths  X  which  are  changing  in  wide 
limits.  Furthermore,  there  along  with  the  attenuation  factor  V  the 
attenuation  factor  for  the  field  of  a  plane  wave  V-|  is  also  given. 

Below  results  are  given  from  the  calculation  of  the  attenuation 
factor  for  arbitrary  values  of  adduced  heights  in  the  area  of  the 
shadow  and  semishadow,  and  the  parameter  q  takes  values  of  q=0  and 
q=  cc  .  Thus  we  encompass,  first  of  all,  the  propagation  of  radio 
waves  having  a  horizontal  polarization  for  any  wavelengths  and  soils; 
secondly,  propagation  of  radio  waves  of  vertical  polarization  for 
wavelengths  and  readily  conducting  soils  (q=0),  and  also  of  short 
waves  and  poorly  conducting  soils  (qsoo  ). 

The  main  calculations  were  made  on  five  members  of  the  series 
(2.6),  (2.7),  (2.9)  and  (2.10).  In  individual  cases  up  to  10-15 
members  of  these  series  were  used,  and  for  construction  of  the  con- 
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nection  graph  up  to  20  members  were  taken.  For  large  and  small  values 
of  adduced  heights  and  distances  different  limiting  formulas  were 
used  ( §  §  3  and  4).  In  the  area  of  the  shadow  the  tables  and  graphs 
are  reduced  to  those  values  of  X,  <  or  Z  and  y  ),  for  which 

in  the  series  (2.6),  (2.7),  (2.9)  and  (2.10)  it  is  permissible  to 
take  only  the  first  term  alone. 

For  calculations  based  on  monomial  formulas  we  present  the  tables 
and  graphs  of  the  factors  entering  into  them,  so  that  calculations 
using  these  formulas  no  longer  present  a  difficulty. 

3.  Attenuation  Factors  in  the  Case  of  Large  Values  of  y 


We  will  consider  the  formulas  for  the  attenuation  factor 
V(  y.  ,yt  ,ji3Jq)  in  the  case  of  large  values  of  one  or  both  cited  heights 

y,  and  Hz  • 

For  example,  under  the  condition  ^>>1  the  expression  for 
V(3r,^;,tj£,  q)  can  be  simplified  considerably,  reducing  it  to  the 
function  only  of  two  variables  and  the  parameter  in  the  following 
manner.  In  the  case  of  large  ’X  and  and  under  the  condition  that 
the  variable  Z  =  T*-  is  not  a  large  negative  number,  the  main  sector 
of  integration  in  (1.8)  corresponds  to  values  of  t  of  an  order  of 
a  unit.  And  since  is  great,  then  to  the  factor  T0(t-yz)  entering 
into  the  function  F(t,j;,^  ,  q)  it  is  possible  to  apply  the  asymptotic 
expression  (see  (8.03)  in  work  [5]) 


't  "t 

r(t—y)»t  Cy  —  0  e 


(3-D 


and  approximately  with  an  accuracy  to  terms  of  an  order  of  -i— 


'7  ~r 


I  i  j’*- I  Vft+t - y— 

(■+*>• r  :  *" 


(3.2) 


or  with  an  accuracy  to  — 

/y 


-r  4  *%-*V7i 

»(/->) -e  y  eT 


(3.2a) 
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Substituting  (3-2a)  into  (1.8),  we  obtain 

't»%  * rz 

V(x,yuyt,q)-t  y  yVifayul). 
where  it  is  assumed  that 

z-x—  )/>v 


(3.3) 


(3.4) 


Thus  the  function  Vi,  being  the  attenuation  factor  for  the  field 
of  a  plane  wave,  simultaneously  makes  it  possible  to  calculate  the 
attenuation  factor  of  a  spherical  wave  q)  for  large  values 

of  ijj,  ,  if  only  Z  =  X is  not  a  large  negative  number. 

Let  us  note  that  function  Vi  in  formula  (3-3),  in  contrast  to 
the  case  when  Vf  enters  into  the  formula  for  the  field  of  a  plane 
incident  wave,  can  also  be  used  for  a  nonuniform  atmosphere  with 
a  linear  course  of  the  refractive  index,  if  in  Z,  and  q  instead 

of  geometric  the  effective  radius  of  the  Earth  a*  is  introduced. 

It  is  more  accurate  to  replace  the  evaluation  ^>>1  of  aPPlica“ 
bility  of  formula  (3-3)  by  the  condition 


(3-5) 


However,  in  the  area  of  deep  shadow,  there  where  the  first  term 
of  series  (2.2)  is  main,  it  is  possible  to  obtain  a  correlation  be¬ 
tween  V  and  Vi  with  a  weaker  limitation  on  32-  In  this  area  it  is 
possible  to  use  the  concept  (3*2)  to  series  (2.2),  preserving  in 
it  terms  of  an  order  of  -J- 

vy 

Designating 


and  j- 


1 

“T 


Mx{y,q)-y  (l  +  -^4)e 

"T  4 «  vjT  +  "iy 
*>y  e 


(3.6) 


we  obtain 


V(x,yx.yvq)' 


,  *  a.*1* 
/-  T 


Ml(y»q)Vl(z.y»<i). 


(3.7) 


In  this  formula  it  is  required  only  that 


<3.8) 


but  then  it  is  applicable  only  in  the  area  of  deep  shadow  (let  us 
say  when  7  -  V7,  >D. 

Formula  (3-3)  is  also  suitable  in  the  area  of  semishadow,  if 
only  Z  is  not  a  large  negative  number  (virtually  when  Z  7-2),  but 
imposes  a  stronger  limitation  on  vj2  . 

If  in  function  V(%  q)  both  and  are  great,  then  it 

is  expressed  approximately  through  function  V 1 i ,  determined  by  the 
formula 


Vu  C.  «*.  1*)+  V»C.q)+  q),  (3.9) 

where 

v<’>c,?)- 


~‘n  i  c 

e  V*\ 


t  M  —  i 

T^(0^e  3qv(0dtt 

-‘T 

w2'(0+e  gwjt) 

(3.10) 


.  - 


~'T  1  “  *•* 

>  M  e  *  <c>°>> 


(3.12) 


'f  -  w.  1  f 
m  7fJe  * 


G<0) 


and 


*-1*. 


.-l/ ^ 


(3.13) 

(3.14) 
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Thus  function  V^,  being  a  function  of  three  variables,  is  ex¬ 
pressed  in  the  form  of  the  sum  of  three  functions,  each  of  which 
depends  only  on  two  variables.  It  is  noteworthy  that  the  term 
does  not  depend  on  the  electric  properties  of  the  soil  and  is  ex¬ 
pressed  through  the  Fresnel  integral,  which  is  used  widely  in  the 
mathematical  theory  of  diffraction.  Below  we  will  examine  the  phys¬ 
ical  essence  of  this  circumstance. 

The  attenuation  factor  V(x,  y  t  q)  is  expressed  through  the 
function  V-|i  in  the  following  manner:  when 


£>o 


V(j:  ,yt,yt.q)- 

*  *|-<*  +  *1  — 

l/xe  (yjj'j)  V,i 


and  when  £  <  0 


f  j rf  !  .  fv.-v.r  1 

Tir *  7 


V{x.yi,y3>  7)  = « 

i  *  i.  *  •  *.*'•»  _  2. 

lyiji  +»  )  -7 

+  l/jce  (ytyj  V'  „  (C.  [*,?). 


(3.15) 


(3.16) 


where  the  variable  £  is  defined  by  the  correlation 


C-x— / Vj-Ky*. 


(3-17) 


and  the  positive  values  of  £  correspond  to  the  shaded,  and  the 

negative  -  to  the  illuminated  areas. 

Analogously  with  large  values  of  reduced  height  y  the  function 

V i  is  expressed  through  the  function  Vn,  and  when  £>0,  in  the 

shadow  we  have  •  , 

ly  f*  -«l 

Vy(x,y,q)^e  y  *VU&\ 


(3-18) 


and  when  £  «<  0,  in  the  illuminated  area 


VM*.y.7) 
1 

T 


ifr*  -4 


+  e  y  Vu  £.*.{)• 


(3.19) 


In  these  formulas  by  T  and  j*.  should  be  understood  the  magnitudes 


(3.20) 


'=*- Vy* 


(3.21) 


Parameter  m  under  the  conditions 


y,  »  i,  yt>  'or’  y  >  1 


(3.22) 


also  should  be  large,  although  for  the  practical  realization  of  this 
values  ij.  and  ft  should  be  very  great.  Thus,  for  example,  if 
20,  then  yU  =  1.5,  and  if  ^  ^  =  1 0 1 ,  then  yut=7.2. 

If  the  following  condition  is  fulfilled 


hi  •!>!>*. 


then  we  have  the  asymptotic  formula 


re“rfs--Sr-  (’<°) 


(3-23) 


(3-24) 


and  then 


—  e 


_1 _ 1_ 

2  /n  C 


(3.25) 


In  this  case  the  dependence  of  function  Vn  on  /U  vanishes 
and  it  becomes  a  function  only  of  the  variables  X  and  q: 


-«  -~J=;+Vw<Z.qnVnC,q). 

*  ,  2  V*  ' 


(3-26) 


Formally  Vi  i  (  £  ,  q )  is  the  value  Vn  (  t  q )  when  f*.~oo 

Vn&q)-VnZ.oo,q)  (3>2?) 

If  f  >  0  (area  of  shadow),  then  the  function  Vn(£>  d)  along 
with  formula  (3.27)  can  be  calculated  based  on  the  series  of  deduc¬ 


tions 


|  **  •  .  *M|  .  * 

V'aG.V)-*  2V*  f  •»(*,) 
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(3.28) 


which  is  obtained  formally  from  formula  (2.3),  if  in  it  instead  of 
H  )  its  asymptotic  expression  (3-2a)  is  taken.  Actually  for 
the  applicability  of  formula  (3-28)  it  is  required  that  the  condition 
(3-23)  take  place. 

Formulas  (3-9)— (3-21 ),  giving  the  solution  to  the  problem  of 
finding  the  diffraction  field  for  large  values  of  reduced  heights, 
were  obtained  byyV.  A.  Fok  [10].  The  conclusion,  for  example,  of 
formulas  (3-18)  and  (3-19)  are  based  on  the  fact  that  in  the  integral 


3 

for  V-j  taken  on  the  broken  path  f"  on  the  sector  from  t=  C  c& 
to  t=0,  for  the  function  $  (t ,  y  ,  q)  the  expression  (1.10)  is  taken, 
and  on  the  sector  from  t=0  to  t=  aO  -  the  expression  (1.11). 

The  first  terms  of  each  of  the  integrals  are  joined  into  one  expres¬ 
sion,  which  is  converted  to  the' form 


«• 

~r 


-  b 


t 


dl 

tx  +  S  ’ 


(3.29) 


where  the  integral  is  taken  based  on  the  infinite  contour  from 

‘  *jr 

t  =  C  *  OO  to  t=o O  ,  enveloping  the  point  5  =-  from  above. 

In  the  case  of  large  values  of  y  this  integral,  calculated  by  the 
method  of  stationary  phase,  gives  that  term  in  the  expression  for 
Vi  which  is  connected  with  the  function  v(°).  In  the  illuminated 
area  this  integral,  furthermore,  gives  the  incident  wave  which  is 
written  out  separately  in  formula  (2.19)  (for  Vi  this  is  a  plane 
wave) . 

The  term  connected  with  the  function  V^D  is  obtained  in  the 
case  of  large  values  of  y  (when  instead  of  ay(t-y)  it  is  possible 
to  take  expression  (3- 2a))  from  the  integral  from  the  second  term 
of  function  <P  on  the  left  half  of  the  circuit  /"  ,  and  the  term 
connected  with  V^2)  «  from  the  integral  from  the  second  term  of 
the  function  £  on  the  right  half  of  the  contour  r  . 

In  a  similar  manner  formulas  (3-15)  and  (3-16)  are  obtained 
for  the  function  V. 

Now  we  will  analyze  the  physical  sense  of  formulas  (3.15), 
(3-16),  (3-18)  and  (3.19).  We  will  begin  with  the  two  latter  ones. 
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With  fulfillment  of  condition  (3.23)  the  term  V^O)  in  expression 
(3-9)  for  V 1 1  is  main  in  the  case  of  finite  values  of  ju  £  ,  since 
in  order  of  magnitude  it  is  yU  times  greater  than  V^*1)  and  v(2). 
Since  ju  increases  very  slowly  with  an  increase  of  y,  and  jjz  »  then 
v(0)  can  be  considered  main  in  comparison  with  V^)+V^2)  only  in 
the  extreme  case.  However,  y(1)+v(2)  change  very  slowly  and  there¬ 
fore  their  sum  can  be  considered  as  a  background,  on  which  the  dif¬ 
fraction  picture  given  by  v(°)  is  superimposed.  Setting  this  back¬ 
ground  off  to  the  side,  we  can  write  the  following  expressions  for 
Vi: 


when  £>0 


v\(z.y,  ?)  =  e 


*<f  >*-*) 


(3-30) 


and  when  £  0 


*»e 


'(t 


—  e 


~'r 


V 


e  ds. 


(3.3D 


where 


(3-32) 


In  the  area  where  these  formulas  are  applicable,  the  field  in 
the  case  of  diffraction  of  a  plane  wave  on  the  spherical  terrestrial 
Earth’s  surface  in  the  point  of  observation  P  (Figure  4)  is  obtained 
with  an  accuracy  to  almost  a  constant  term  the  same  as  in  the  case 
of  diffraction  on  a  screen  passing  through  the  center  of  the  ter¬ 
restrial  sphere,  and  the  edge  of  this  screen  is  the  line  of  inter¬ 
section  of  the  plane  of  the  horizon  for  point  P  in  the  geometric 
boundary  of  the  shadow  of  the  incident  plane  wave  (Figure  4). 

Actually  we  will  designate  through  the  angle  of  diffraction, 
i.e.,  the  angle  between  the  direction  of  propagation  of  the  incident 
plane  wave  and  the  direction  of  the  radius-vector  r=TP,  drawn  from 
the  edge  of  the  screen  perpendicular  to  it  in  the  point  of  observa¬ 
tion  P;  in  this  case  we  consider  'lf,>  0,  if  P  lies  in  the  shaded 
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area  (as  in  Figure  4),  and  in  the  opposite  case  ^  0.  Then  based 

on  the  theory  of  diffraction  on  a  half-plane  the  fields  are  expressed 
through  the  function 


~lT  - 


y/T  J 


(3.33) 


~  ‘  4  1 


/*r  Ikr  coti  e 

U  —  e  — e  — 


p  *  ight 

■=—  l  e  is  (cbct), 


(3.34) 


where 


t  \j2kr  sin  . 


(3-35) 


u  > 

n  ^ 


Figure  4.  Equivalent  screen  for  a  plane  wave 


Simple  geometric  considerations  on  the  basis  of  Figure  4  lead 
to  the  fact  that  for  small  angles  of  diffraction  ^ ,  large  values 
of  ij  and  finite  £  ,  we  have  the  correlations 

kr  rr,  2/Vf’  \/7 

-vy-.  I  (3-36> 


Thus  the  Fresnel  integrals,  through  which  the  attenuation  factor 
in  formulas  (3.30),  (3-31)  is  expressed,  have  the  same  arguments 
as  the  diffraction  field  of  the  screen  in  Figure  4. 


The  factor  e~c*i  in  front  of  the  integrals  in  (3-30)  and  (3.31 ) 

'L 

is  equal,  as  is  readily  checked  from  (3*36),  to  the  factor  e"  2 
obtained  in  the  case  of  using  the  approximate  formula 

«$>* 

cos  -V  =  1 - — 

The  term  e  'T'  *  ,  giving  the  incident  plane  wave  in  formula 

(3-31),  corresponds  to  the  term  e in  formula  (3-3*0.  The 
remaining  differences  in  formulas  (3*30),  (3^30  and  (3-33),  (3-3*0 
are  connected  simply  with  a  different  reading  of  phases  in  functions 
Vi  and  U. 


Figure  5.  Equivalent  screen  for  a  spherical  wave. 

The  same  physical  essence  is  held  by  the  term  jn  the  ex¬ 

pression  for  the  attenuation  factor  V  [formulas  (3-15),  (3-16)]. 

In  this  case  the  distances  TQ=ri  and  TP=r2  (Figure  5)  are  connected 
with  the  cited  heights  of  the  corresponding  points  by  the  correla¬ 
tions 


*Tj  *■  2Af  *  '\f~Ti- 


(3-37) 


The  screen  itself  is  defined  as  a  half-plane,  passing  through 
the  center  of  the  Earth  and  having  as  its  own  edge  the  straight  line 
of  intersection  of  the  horizontal  planes  for  the  corresponding  points 
P  and  Q.  The  angle  of  diffraction  *P  is  defined  as  the  angle  between 
these  planes  and  is  considered  positive,  if  P  is  not  found  in  the 
zone  of  direct  visibility  of  point  Q.  Correlations  (3-37)  are  valid 
under  the  assumption  that  angle  V*  is  small. 
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As  is  known  from  the  classic  theory  of  diffraction  the  field 
from  a  nontransparent  screen  of  such  a  form  is  expressed  (at  small 
angles  of  diffraction,  which  only  interest  us)  through  the  Fresnel 
integral  from  the  argument 

.  /  2/rr,r,  •!» 

\  tttt:  2 '  (3.38) 


Since  as  before  the  angle  of  diffraction  is  connected  with 
the  variable  X  by  the  correlation 


t  =  M  *, 


then  formula  (3-38)  gives 


»  *  y\  +  yyj 


(3.39) 

( 3  -  40 ) 


which  coincides  with  (3-13)  and  (3-14). 

Thus  the  term  in  the  expression  for  the  function  V 1 1  ( S' ,  q 

corresponds  to  the  standard  Fresnel  diffraction,  which  according 
to  the  Huygens  principle  is  obtained  from  the  nontransparent  Earth 
the  same  as  from  a  nontransparent  screen,  positioned  so  that  it  pro¬ 
duces  (for  the  given  positioning  of  the  corresponding  points)  the 
same  geometric  boundary  of  shadow  as  the  Earth. 

The  diffraction  field  is  obtained  according  to  the  Huygens  prin¬ 
ciple  as  a  result  of  the  effect  of  that  part  of  the  incident  plane 
or  spherical  wave  which  for  the  given  point  of  observation  P  is  not 
obstructed  by  a  nontransparent  obstacle. 

In  previous  works  on  the  theory  of  diffraction  propagation  of 
radio  waves  around  the  surface  of  the  Earth  repeated  attempts  were 
made  to  apply  the  Huygens  principle  and  to  establish  a  connection 
between  diffraction  around  the  Earth  and  diffraction  on  the  edge 
of  a  screen.  However,  those  elementary  considerations,  which  were 
expressed  in  this  connection,  did  not  make  it  possible  to  determine 
the  boundaries  of  applicability  both  of  the  Huygens  principle  and 
of  those  calculation  formulas  of  the  type  (3.33),  ( 3 • 34 )  which  are 
obtained  with  its  help  and  coinoide  with  formulas  of  diffraction 
on  a  wedge  or  half-plane. 
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In  particular  these  formulas  were  applied  to  the  case  of  a  deep  shadow 
(large  positive  X  )>  where  they,  according  to  formula  (3.25),  give 
an  attenuation  of  the  field  inversely  proportional  to  £  >  while 
according  to  strict  theory  (formula  3.28)  the  field  fades  with  an 
increase  of  £  exponentially.  Therefore  it  is  not  surprising  that 
the  Huygens  principle  gave  a  field  which  is  many  times  stronger  than 
that  actually  observed.  Even  in  the  comparatively  recently  published 
book  by  Bremmer  [9]  in  the  interpretation  of ‘this  interesting  and 
important  question  total  helplessness  is  revealed.  Only  in  the  works 
by  V.  A.  Fok,  stemming  from  the  boundary  integral  as  the  main  expres¬ 
sion  for  the  attenuation  factor  V  (and  not  from  a  series  of  deductions 
as  this  is  usually  accepted),  this  question  was  cleared  up,  which 
then  made  it  possible  to  effectively  calculate  the  attenuation  factor 
for  high  values  of  cited  heights  y  [10]. 

From  a  comparison  of  formulas  (3-26)  and  (3-28)  it  is  evident 
that  only  in  the  deep  shadow  the  term  V^°)  is  not  main;  v^^+V^2) 
has  the  same  order  of  magnitude,  and  the  sum  of  all  three  terms  has 
the  form  (3.28),  i.e.,  with  it  decreases  considerably  more 

rapidly . 

In  general,  as  we  already  mentioned,  the  term  v(°)  in  the  expres¬ 
sion  for  V i  •(  is  main  only  for  high  values  of  y  .  Even  when  y  =1000 
m-JT  =5-6,  so  that  when  X  =0  the  term  V^0)  is  only  10  times 
greater  in  absolute  magnitude  than  the  summ  V^^+V^2).  Therefore 
it  can  almost  never  be  disregarded.  The  term  V^°)  can  be  considered 
main  only  because  at  low  values  of  X  it  gives  a  rapid  course  of 
the  functions  V-|i,  V-j  and  V:  in  the  case  of  positive  values  of  X 
thanks  to  this  term  the  functions  V-|  and  V  increase  rapidly  from 
low  to  finite  values,  and  with  negative  values  of  J  -  they  oscillate 
rapidly . 

The  terms  and  V^2),  depending  on  the  electric  properties 

of  the  soil,  determine  the  influence  of  the  material  and  form  of 
the  "nontransparent  obstacle"  -  the  Earth,  on  the  diffraction  field 
-  the  circumstances,  which  the  classic  theory  of  diffraction 
( Fresnel -Kirchhof f )  does  not  take  into  account.  In  addition  to  the 
parameter  q  these  terms  depend  only  on  J  and  therefore  at  low  values 
of  X  can  be  considered  almost  constant. 
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Let  us  note  that  the  function  V-|i(J,  q)  [see  (3*26)3  emerges 
in  the  solution  of  problems  of  diffraction  not  only  on  a  sphere, 
but  also  on  bodies  of  another  form.  Thus  these  functions  were  used 
in  the  work  of  A.  S.  Goryainov  [11],  who  considered  diffraction  on 
an  infinite  round  cylinder.  In  this  work  the  following  designations 
are  used : 


fO  -  -  Vu  C,  oc) 

7(0 --[IW&00)+ 00)1 
7(C)-  -  [V'<"(C.O)+  V' <=>(;,  0)1 


4.  Attenuation  Factors  at  Low  Values  of  y 


In  the  previous  section  we  considered  formulas  for  the  calcula¬ 
tion  of  attenuation  factors  in  the  case  of  large  values  of  y  .  The 
case  of  low  values  of  y  should  also  be  considered  in  particular 
and  the  formulas  necessary  for  calculations  given.  The  difficulty 
here  is  that  at  low  y  the  area  of  half-shadow  has  corresponding 
low  values  of  X  and  2  ,  at  which  the  series  of  deductions  converge 
very  slowly. 

Therefore  here  it  is  necessary  either  to  calculate  the  atten¬ 
uation  factors  directly  based  on  their  integral  presentations,  using 
for  simplification  of  the  integrands  the  fact  that  y  4^1 ,  or  to  cal¬ 
culate  in  series  of  deductions  the  remainder  terms,  which  in  the 
area  of  half-shadow  at  small  y  no  longer  can  be  disregarded. 

Let  us  note  first  of  all  that  if  function  )  satisfies  the 
equation 


r(y)-v-y)fiy) 

with  the  initial  conditions 


(4.1) 


/(0)  =  1.  /’(O)  “  — 


26. 


(4.2) 


I  ... 


then 

f(y)  **  ch  (yV~t) - -  sh(y/7), 

/  <  (.4.3; 

if  only 

^ri<1  <*-'•> 

(see  [5],  pages  45-46).  y) 

It  is  easy  to  see  that  the  height  factor  ir(ft)  , 

entering  into  the  series  of  deductions  (2.2)  and  (2.3),  satisfies 
equation  (4.1)  with  the  initial  conditions  (4.2). 

Therefore 


w(0  Vi, 


(4.5) 


under  the  condition 


2|r7;} 


<  1. 


(4.6) 


If  and  satisfy  the  condition  (4.6),  then,  substituting 
(4.5)  for  and  y2  in  (2.2),  we  obtain 


'T 


V[x,y„  y»,  «»)*-  — e  2  V  X 

*  ut  *  _ 

x  —  sh Cr, /l?)  sh (yr » '  *;')  (4-') 

••1 


V(*,.yi,y2.0)  =  e  2  V~x  .< 


Ixt' 


X  v  -e  ch  (y,  /  l/)  ell  (y2  /  ).  (4.8) 


l»l 


Introducing  the  designations 


♦  (•».’!)  =  «  ,  i - ? 


•  •I 


t—  -  "  irt/tWV 

'^E-Sr1 


••1 


(4.9) 


(4  10) 
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*U*e  12 - 


2 7r  * 

r-r  3/,;  .  5x_i»  “'‘S' 


(4.15) 


— >»W  - 


,  '» 

“2/re  * 


(4.16) 


If  at  low  values  of  and  <jz  the  value  of  %  is  not  small 
(<■ V^.1),  then, since  the  factor  when  *Vx  in  (4.15)  and  (4.16)  is 
a  number  of  order  3,  the  integrals  in  (4.13)  and  (4.14)  are  a  small 
correction  to  the  first  N-1  members  of  the  series  and  they  can  be 
disregarded.  In  the  opposite  case  these  integrals  have  the  same 
order  as  the  sums. 

From  such  considerations,  as  was  pointed  out  by  V.  A.  Fok, 
for  small  values  of  4j  it  is  possible  to  obtain  the  formula 

Vt  (zj.q)  - 

,  fwch(/0-iili(^) 

.  ‘  e  - — - dt, 

V*r  vf(t)  —  qw(t) 


or,  dividing  the  contour  f  into  parts,  as  in  §3, 

V,(aj,9)  = 


ch(j>/7)~  -y  si  Kyi  t  ) 
w'  (t)  —  qii'U) 


-  ..gt  J  <>cli(K)— e  J  ~-shO  • 


+  4=  c 


w/(0+e  q*\U, 


(4.17) 


(4.18) 


where 


r=y/7e'  3. 


Using  the  last  formula  calculations  were  made  for  graphs  of 
connection  in  Figures  6-9.  However,  this  formula  gives  approximately 
the  same  result  as  the  much  simpler  expressions  (4.23),  (4.24), 
which  are  given  below. 

We  will  also  use  the  following  formulas • 

A.  y,  <  l 

-  V  (x,yu  y,,  oo)  , yx  <L¥(x,0,  y,.  oo) 

dy{ 

v  (x>  y>.  y».  0)  -  V  (x,  0.  0).  (4.20) 

B.  yi  <  1,  y2  «  1 

V  {x,y „y2,  oo)-y,  y..  2LV!*’.  °»  ^°>(4  21) 

rfy,  dys 


V  {. x,yuyt,  0)  «  V  (x,  0. 0. 0)  (4.22) 


.  C.  y  «  1 


V,  (z,  y,  «,) -j,fL£fr°»  °°>  «y  /(_-) 

(4.23) 

V'l  (X, y,  0)  -  V, (z,  0.0)  =  * (z).  (4.24) 

Me  where 

<?^(z,0,.y.,  y)  _ 

<?yi 


-i* 


4  n  /■  ~~  _ _  ix  t  »* 

e  S/wJe  *,(y.),  (4.25) 


«-i 


<?*  V(x,0,0,oo)  _  —  2^  mV 

- -  ;  1  -e  2  >  V  e  (4-2H) 

vi  ®ys  ,tr, 

dy 


•  i*tt'  .  c  «*• 

i  2/«S- - -4*1— — dt  (4J27) 

,-i  »'((/•)  /«  t  w(() 

^(z)»K,(z,0.0)  = 

-i2»'«2  — - -»“4=4  ~~T~dt.  (4.2S. 

V«i  w(t) 
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(Functions  £  (1)  and  g(Z)  play  an  important  role  in  a  number 
of  other  problems  and  have  been  tabulated  in  detail. 

Let  us  note  that  formulas  (4.21)  and  (4.22)  can  be  obtained 
also  from  formulas  (4.7)  and  (4.8),  if  in  all  essential  terms  of 
series  (4.7)  and  (4.8)  it  is  permissible  to  replace  the  hyperbolic 
functions  by  the  first  terms  of  their  expansions  into  a  power  series. 
An  analysis  of  this  requirement  gives  the  following  condition  of 
applicability  of  formulas  (4.21)  and  (4.22)’ 


Vx 


<  1. 


(4.29) 


In  a  similar  way  formulas  (4.23)  and  (4.24)  can  be  derived  from 
(4.18),  and  the  condition  of  their  applicability,  obtained  from  the 
requirement  of  the  possibility  of  replacement  on  the  main  sector 
of  integration  of  the  hyperbolic  functions  by  the  first  terms  of 
their  expansion,  has  the  form 


>1*1  <1  — *»1  \  (4.30) 

y  C  1  I*!<1  *  *>0  (’ 

fcft. 

5.  Reflection  Formulas.  Joining  with  Diffraction  Formulas 

As  is  known,  for  the  calculation  of  field  intensity  in  the  illum¬ 
inated  area  at  relatively  small  distances  from  the  source  reflec¬ 
tion  formulas  are  used.  These  are  based  on  the  concept  that  the 
reflection  of  beams  from  the  flat  surface  of  the  Earth  according 
to  the  laws  of  geometric  optics  is  such  that  an  electromagnetic  field 
in  a  given  point  of  the  illuminated  area  emerges  as  a  result  of  the 
interference  of  the  "direct”  and  "reflected"  radio  beam.  As  theory 
shows,  for  the  ideally  reflecting  flat  Earth  this  reflecting  treat¬ 
ment  gives  completely  accurate  results  (with  any  polarization  of 
the  propagating  radio  waves),  at  the  same  time  that  for  the  flat 
Earth,  possessing  a  finite  conductivity,  the  reflecting  formulas 
are  only  approximately  correct,  and  in  the  case  of  glancing  incidence 
of  the  beams  for  waves,  of  vertical  polarization  for  example,  should 
be  replaced  by  Weil-Van  der  Pol  formulas. 
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The  application  of  the  stationary  phase  method  to  contour  in¬ 
tegrals  (1.8)  and  (1.17)  for  V(X,  q)  and  V-|(z,<f,  q)  makes  it 
possible  to  obtain  in  the  illuminated  area  for  these  functions  the 
following  formulas  (see  [1]  and  [5]): 


—  in  -T-  *“* 


(5.1) 


(5.2) 


These  are  reflection  formulas  for  the  attenuation  factors  V 
and  Vi,  taking  into  account  the  curvature  of  the  Earth  and  genera¬ 
lizing  the  known  reflection  formulas  for  a  flat  Earth.  The  first 
term  in  formulas  (5.1)  and  (5.2)  gives  the  incident  wave  (spherical 
for  V  and  plane  for  V-|),  the  second  term  -  reflected  wave:  co  or 


LL 

Y  is  its  phase,  % ~  coefficient  of  Fresnel  reflection 
[for  glancing  incidence,  since  only  for  this  case  are  formulas 
(1.1)— (1.4)  and  ( 1 . 1 6 )  valid],  -/T  and  V*  -  factors,  taking 
into  account  the  additional  convergence  of  rays  as  a  result  of  their 
reflection  from  the  spherical  surface  of  the  Earth. 

The  parameter  p,  designated  in  like  manner  in  (5.1)  and  (5.2), 
and  having  the  same  geometric  sense  in  both  cases  (which  we  will 
speak  of  below),  is  determined  for  V  from  the  equation 


/  y,  +  v 'y3H  P'  *  2 p  -i  x, 


and  for  V-|  according  to  the  formula 


p~±(r'z>  +  3y-2z). 

w 


(5.3) 


(5.4) 


The  solution  of  equation  (5.3)  can  be  presented  in  the  form 


x  +  z 


where 


z=2?sin  —  , 


(5.5) 


(5.6) 
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t5 . 7 ) 


p-  =  —  (jc‘  +  2y ,  -|  2y,)  (p  >  0), 
3 


-v (3*1  —  >'-•)  (  ~  -■ ■  „  •  -  jL\ 

s,m* 

(5.8) 

Further 

t  4x  2  ,  12 

«*•  =  —  3 p*x+  2/7  (y,  +y,  -  x*)  + 

(5.9) 

+  x(y1+yi)~  j*'. 

(5. 10) 

A_,  P*  . 

3/ix  +  jr»—  yt-y. 

(5.11 ) 

For  Vi  we  have 

the  correlations 

*• 

(5.12) 

(5.13) 

—  Vr  +  Zy.  . 

(5.14) 

The  geometric  sense  of  the  parameter  p  lies  in  the 
in  characterizes  the  area  in  which  the  point  j',  or 

found . 

If  T  is  the  angle  of  incidence  of  the  beam,  then 

fact  that 
Zfy  is 

p  —  M  cos  i. 

(5.15) 

from  which  it  is  immediately  clear  that  p=0  if  the  point  lies  on 
the  geometric  boundary  of  the  shadow,  p>0  if  the  point  lies  in  the 
area  of  light,  and,  finally,  p^O  if  it  lies  in  the  shadow.  Let 

us  stress  that  this  important  factor  does  not  depend  on  the  properties 
of  the  soil. 

For  V •)  i  it  is  also  possible  to  obtain  a  "reflection"  formula, 
i.e.,  a  formula  which  is  applicable  in  the  case  of  large  negative 
t  •  For  this  in  the  second  member  of  formula  (5.2)  ( V i i  in  the 
area  of  light  by  definition  does  not  contain  an  incident  wave) 
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it  is  necessary  to  introduce  in  place  of  z  the  variable 

and  to  consider  that  1  and  y  are  great,  and  £  is  finite  (and 

negative ) . 

Then  the  result  is  that  Vn  does  not  depend  on  ^  and  has  the 


form 


q  +  ip  2 


12  ’ 


(5.16) 


where 


(5.17) 


Let  us  note  that  the  reflection  formulas  for  V,  Vi  and  V-ji  con¬ 
vert  one  another  in  the  case  of  large  values  of  ij  in  the  same  man¬ 
ner  as  the  diffraction  formulas  for  these  functions. 

Reflection  formulas  of  the  type  (5.1)  frequenly  are  used  for 
calculations,  considering  in  them  Va  si,  i.e.,  disregarding  the  addi¬ 
tional  divergence  of  beams  as  a  result  of  their  reflection  from 
the  Earth. 

Figure  6  gives  an  example  of  curves,  calculated  using  formula 
(5.1)  (solid  curves)  and  in  the  case  of  V"a=1  (broken  curves).  As 
it  should  be  expected,  the  latter  oscillate  more  sharply  and,  in 
contrast  to  the  first,  reach  in  the  minimum  of  zero  the  same  position 
of  the  maximum  and  minimums  for  both  curves. 

As  theoretical  investigations  show,  the  reflection  formulas 
are  applicable  in  the  illuminated  area  under  the  condition  that 
the  parameter  p  is  positive  and  great. 

For  refinement  of  how  great  the  parameter  p  should  be  in  order 
that  the  reflections  would  give  a  specific  accuracy,  we  made  a  series 
of  calculations  on  the  joining  of  reflection  formulas  with  diffrac¬ 
tion  . 

Figures  2-5  give  the  results  of  investigations  of  functions 
» fh ,e>^)  and  by  diffraction  series  (2.6)  and  (2.9) 

and  by  refelction  formula  (5.1)  for  =0 . 4 ,  0.8,  2.4,  4.8. 

The  oscillating  parts  of  the  curves  are  constructed  according 
to  the  reflection  formula.  The  values  of  parameter  p  are  plotted 
on  them>.  The  figures  in  brackets  show  on  what  number  of  terms  of 
series  (2.6)  or  (2.9)  the  given  points  are  calculated.  The  little 
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circles  denote  the  points  which  are  calculated  based  on  more  than 
five  members  of  the  series. 

In  analyzing  these  drawings  we  see  that  already  with  p=0.7  for 
q  =  oO  ,  and  for  q=0  even  when  p=0.5  there  is  an  approximate  joining 
of  curves  constructed  by  reflection  and  diffraction  formulas. 
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Figure  6.  Influence  of  additional  divergence  on  curves  constructed 
using  the  reflection  formula. 


This  result  shows  that  reflection  formulas  actually  have  a  wider 
area  of  applicability  than  this  would  be  expected  on  the  basis  of 
general  analytical  evaluations. 

In  this  case  for  q=oO  the  joining  takes  place  after  the  passing 
of  the  last  maximum  of  absolute  value  of  function  V,  i.e.,  in  the 


area  of  its  monotony,  where  it  is  possible  yet  to  be  limited  to  five 
members  of  the  series  of  deductions  (2.6)  (for  not  too  small  X  )• 

For  q=0  (also  for  not  too  small  ^  )  the  joining  takes  place 
on  the  last  maximum  of  function  V  when  p<*  0.5-0. 7  and  also  in  the 
area  of  suitability  of  the  five-term  formula.  However,  the  circum¬ 
stance  that  the  joining  takes  place  on  the  maximum  hampers  the  graphic 
representation  of  the  results  of  calculations  based  on  a  five-member 
series  (see  §6  below).  The  joining  of  reflection  and  diffraction 
formulas  for  the  attenuation  factor  Vi  is  given  in  drawings  6-9- 

Here  the  series  of  deductions  (2.7)  and  (2.10)  do  not  make  it 
possible  to  obtain  the  joining  with  reflection  formulas  (even  if 
25  terms  are  taken  in  them)  and  we  had  to  use  other  methods.  Thus 
for  y  =4.8  formulas  ( 3 . 1 8 ) -( 3 . 1 9  )  already  turned  out  to  be  suitable 
for  large  values  of  y  .  For  y=0.4  and  ,y=0.01  we  used  formulas 

(4.23)  and  (4.24),  and  for  <y=0.4  also  formula  (4.18). 

In  draw.  6-9  the  points  obtained  with  formulas  (3-18)  and  (3-19) 
are  noted  by  little  circles,  points  calculated  with  formulas  (4.23) 
and  (4.24)  are  noted  by  triangles,  and,  finally,  the  results  of  cal¬ 
culations  by  (4.18)  -  by  small  squares. 

The  broken  curve  in  drawing  8  represents 

V,(z,  0,  0)-g{z), 

-  the  extreme  curve,  to  which  Vi  (*,</,  0)  strives  when  jf-+Q  [formula 

(4.24) ],  although  the  striving  to  the  limit  is  nonuniform  and  depends 
on  the  value  of  2.  Thus  in  the  area  of  shadow  and  semishadow  already 
for  ^  =0.4  the  function  Vi(s,y,  0)  is  very  close  to  g(a),  in  the 

area  of  light  Vi(Z,y,  0)  for  t/=0.4  has  a  completely  different  course 
than  for  g(x) .  However,  the  smaller  the  value  of  y  ,  the  further 
into  the  area  of  light  is  formula  (4.24)  suitable.  Thus,  when 
^=0.01  for  2<  -1.2  the  points,  calculated  by  the  reflection  formula 
(noted  by  crosses),  actually  coincide  with  g(a).  The  same  is  true 
when  q=oo  [formula  (4.23)]. 

These  circumstances  agree  with  the  evaluation  of  (4.30)  given 
in  §4  for  the  applicability  of  formulas  (4.23)  and  (4.24).  As  was 
already  stated  above,  formula  (4.18)  gives  the  same  results. 

The  joining  of  the  reflection  and  diffraction  formulas  for 
V i  (Z,  ,  q)  takes  place  at  those  same  values  of  p  as  for  V(y,  y  t ,  ^t,q ) , 
namely  when  p=0.7  for  q  =  e©  and  when  p-^0.5-0.7  for  q=0.  In  both 
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cases  the  joining  takes  place  on  the  maximum  of  the  absolute  magni¬ 
tude  of  the  function  or  even*  further  to  the  left.  The  joining  of 
Vii  with  formula  (5.16)  is  shown  in  drawings  36-39.  It  takes  place 
when  p  -^0 .7-1 .0  . 

Let  us  note  that  the  joining  of  reflection  and  diffraction 
formulas  for  q)  takes  place  at  positive  or  small  negative 

values  of  ZL  ,  i.e.,  in  the  area  where  formula  (3-3)  of  connection 
between  V  and  Vi  is  applicable.  This  makes’ it  possible  for  us  to 
draw  the  conclusion  that  the  joining  of  diffraction  and  reflection 
formulas  for  V  ( *,  ,  q)  for  the  case  when  y2:s>1  and  y,  is  finite 

or  small  takes  place  also  when  p~0.7.  The  joining  of  Vii  with 
formula  (5.16)  shown  in  drawings  36-39  ma"kes  it  possible  for  us  to 
make  such  a  conclusion  for  ,  »1 ,  ^  »  1 . 

For  determination  of  the  area  of  applicability  of  reflection 
formulas,  i.e.,  the  area  where  p)^0.7,  we  present  special  graphs 
(drawings  10  and  23  -  see  §6). 

6.  Description  of  Graphs  and  Tables 

As  was  already  stated  above,  we  will  give  the  tables  and  graphs 
of  attenuation  factors  V  and  Vi  for  the  values  q»  ao  and  q=0. 

The  calculations  of  the  function  V  for  these  values  of  q  were 
made  for  several  fixed  ratios  of  cited  heights. 

m  „  ll.  -  l;  0,75;  0,5;  0,25, 

y* 

and  values  of  'X  from  0.5  to  6  (q=oe)or  to  7  (q=0)  and  from  0.1 
to  4.8,  and  calculations  were  made  for  all  pairs  (  for  which 

by  using  a  five-member  series  it  was  possible  to  obtain  results 
if  only  with  graphic  accuracy. 

Values  of  y  >4.8  were  not  taken  in  these  calculations,  since 
for  them  it  is  possible  to  obtain  the  attenuation  factor  V  with  the 
help  of  function  Vi  USing  formulas  (3-3)  and  (3-7)  or  (3. 15)-(3- 16). 

For  values  of  X>6  the  functions  V  and  Vi  did  not  have  to  be 
calculated,  since  for  such  T  and  ‘J,  .8  already  the  use  of  only 

the  first  member  of  the  series  of  deductions  is  permissible  (see 
below ) . 
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Results  of  the  calculation  of  function  V  using  the  five-member 
series  are  given  in  tables  1-16. 

The  function  of  attenuation  Vi  was  calculated  also  for  t^=0.1,. 

4.8  and  values  2  from  0.5  to  6  (q  =  ek?  )  or  to  7  (q=0)  (tables 
17-20) . 

Tables  21-24  give  the  first  five  height  factors  and  their  de¬ 
rivatives  for  y=0.1,  0.2,....,  5.0. 

Furthermore  tables  25-33  give  auxiliary  functions  for  calcula¬ 
tions  using  formulas  (4.11),  (4.12),  (4. 19)— (4.24) . 

In  the  tables  given  the  error  can  reach  one  and  a  half  units 
of  the  last  sign  written  out. 

Small  type  is  used  for  the  figures  in  which  tne  error  can  reach 
three  units. 

In  practical  calculations  it  is  most  convenient  to  use  the  graphs 
which  are  appended  to  this  work. 

Drawing  1  gives  the  standard  distance  (A)  and  the  standard 
height  h0{^)  for  obtaining  the  reduced  distances  and  heights  by  using 
formulas  (1.5),  (1.18),  (1.19).  On  the  axis  of  abscissas  in  a  loga¬ 
rithmic  scale  is  plotted  the  length  of  the  wave  in  meters.  On 
the  axis  of  ordinates,  also  in  a  logarithmic  scale,  are  plotted  the 
variables  s0(^)  and  h$(A)  in  meters.  Curves  are  given  for  different 
values  of  effective  radius  a!*  from  Ci-d  to  £*  =  4 a. 

The  basis  for  the  graphs  of  the  attenuation  factor  V  for  q=oO 

(drawings  11-14)  and  q=0  (drawings  15-18)  are  tables  1-16,  in  which 

y 

this  factor  is  given  for  fixed  values  of  the  ratio  m  =  ~j- . 

•'a 

Furthermore  for  these  graphs  all  the  auxiliary  formulas  for  large 
and  small  values  of  y  were  used  and  a  series  of  precalculations 
was  made  for  values  of  m  which  are  not  found  in  tables  (1-16). 

In  the  most  simplest  manner  the  function  V  depends  on  X  in 
such  a  way  that  based  on  this  variable  the  linear  interpolation 
is  permissible.  Therefore  all  the  graphs  were  constructed  for  fixed 
values  of  %  ,  taken  through  0.5. 

On  the  axis  of  abscissas  of  all  these  graphs  is  plotted  the 
variable 
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On  the  axis  of  ordinates  in  drawings  11,  12,  15  and  16  the 
lg  I VI  +  10  or  arc  V  are  plotted.  In  each  sketch  a  family  of  curves 
is  given  for  fixed  values  of  fa.  For  each  of  the  values  of  q  =  o* 
and  q=0  a  drawing  is  given  which  consists  of  two  bands,  on  which 
such  graphs  are  grouped  for  fixed  values  of  which  are  arranged 
in  increasing  order. 

For  convenience  of  interpolation  the  graph  for  *X= 3.5,  with 
which  the  first  band  always  ends,  is  repeated  on  the  second.  The 
broken  lines  which  intersect  the  families  of  curves  represent  the 
line  of  the  horizon  p=0  (or,  which  is  the  same,  ^  —  -/*  =o). 

In  drawings  13,  14,  17,  18  on  the  axis  of  ordinates  ^(less 
than  g  )  is  plotted,  and  on  the  diagram  the  lines  of  equal  values 
of  the  function  are  drawn.  The  broken  curve  gives  the  line  of 
the  horizon,  and  the  dot-and-dash  curve  represents  the  line  p=0.7, 
above  which  the  reflection  formula  is  applicable. 

As  was  noted  in  the  previous  section,  for  q=  oo  the  closing 
of  the  reflection  and  diffraction  formulas  takes  place  to  the  right 
of  the  maximum,  closest  to  the  area  of  the  3hadow,  of  the  modulus 
of  the  attenuation  factor  in  the  area  of  monotonicity  of  the  func¬ 
tion  (see  Fig.  2),  whereas  for  q=0  this  closing  takes  place  on  the 
maximum  itself  or  even  to  the  left.  This  explains  why  for  q=0  the 
graphs  do  not  reach  the  dot-and-dash  line  p=0.7  everywhere,  i.e., 
the  reflection  formulas,  because  the  graphic  representation  of  a 
family  of  non-monotone  functions  is  very  difficult. 

Drawings  19-22  give  the  attenuation  factor  V  for  small  values 
of  %,$ |  and  fa.  For  each  of  the  values  #=0.007,  0.01,  0.03,  0.05, 
0.1  and  0.25  a  separate  graph  is  given,  and  all  these  graphs  are 
assembled  on  one  sheet. 

On  the  axis  of  abscissas  Jfc^O.1  is  plotted,  and  on  the  axis 
of  ordinates  -  the  function  I V I  or  arc  V  and  for  each  value  of 
Uj<0.1  a  separate  curve  is  given.  For  intermediate  values  of  fa 
the  function  can  be  obtained  either  by  interpolation  on  the  graph 
or  by  formulas  (4.11),  (4.12)  with  the  help  of  tables  25,  26. 

For  the  attenuation  factor  Vj(*,^,q)  we  give  graphs  of  three 
types  (for  each  of  two  values  of  q):  dependence  of  lg  I  Vi  I  and 
afrc  Vi  on  2  with  fixed  values  of  y  (drawings  24,  25,  30,  3D, 
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the  dependence  of  these  functions  on  y  with  fixed  values  of  "Z. 
(drawings  26,  27,  32,  33),  and,  finally,  their  level  lines  (drawings 
28,  29,  34,  35).  On  graphs  of  the  last  type  the  variable  Z  is 
plotted  on  the  axis  of  abscissas,  and  on  the  axis  of  ordinates  - 
the  variable  y  .  The  broken  line  (on  all  the  graphs)  denotes  the 
line  of  the  horizon,  the  dot-and-dash  -  the  line  ps0.7  -  the  boundary 
of  applicability  of  reflection  formulas. 

For  Vi,  as  was  already  pointed  out  above,  the  series  of  deduc¬ 
tions  converge  more  poorly,  and  it  is  not  possible  to  achieve 
closing/joining  with  negative  formulas  on  a  mass  scale. 

However,  for  large  or  small  values  of  y  it  is  always  possible 
to  reach  the  areas  P'5.0.7  with  the  help  of  formulas  (3* 1 8 ) — ( 3 . 19) 
and  (4.23),  (4.24).  The  enumerated  graphs  of  the  attenuation  fac¬ 
tors  V  and  Vi  we  preface  with  diagrams  (drawings  10  and  23)  which 
give  the  areas  of  applicability  of  the  reflection  formulas,  i.e., 
the  areas  where  p$>0.7,  and  these  diagrams  are  given  in  reduced 
and  geometric  coordinates. 

In  drawing  10  the  graph  of  applicability  of  reflection  formula 
(5.1)  for  function  V  is  given. 

On  the  axis  of  abscissas  of  this  graph  the  ratio  ***  =  is 
plotted,  and  on  the  axis  of  ordinates  -  yt  .  For  each  value  of  ^ 
(through  0.5)  curves  are  given  on  which  the  condition  p=0.7  takes 
place.  Each  curve  divides  the  area  of  all  (*/,»&)  into  two,  in 
one  of  which  (upper)  at  a  given  value  of  Z  the  reflection  formula 
is  applicable. 

Let  us  note  that  increases  upwards,  and  </t  increases  with 

y 

movement  to  the  left,  when  \m  =  decreases.  Thus  the  upoer  area 

/i 

corresponds  to  greater  heights,  and  the  lower  -  to  lesser. 

Since  the  distance  and  height  are  connected  with  the  reduced 
coordinates  by  the  correlations 

*-«•(■«•)<* 

then,  in  order  to  switch  from  ne  and  y  to  $  and  k  ,  it  is  sufficient 
to  change  the  scale  on  the  axis  of  ordinates  and  to  replace  the  values 
of  X  on  the  curves  in  accordance  with  the  indicated  formulas.  On 
our  graph  scales  of  geometric  heights  (in  meters)  and  distances  (in 

40. 


kilometers)  are  given  for  several  wavelengths  A  and  effective  radius 
of  the  Earth  £*=10,000  km. 

Drawing  23  gives  the  areas  of  applicability  of  reflection  formu¬ 
la  (5.2)  for  function  Vi-  Here  parameter  p  depends  only  on  two  var¬ 
iables  Z~  and  y  .On  the  axis  of  abscissas  the  variable  2.  is  plotted, 
and  the  distances  (in  kilometers),  corresponding  to  several  wave¬ 
lengths  "X  ,  of  the  observation  point  from  the  boundary  of  the  geome¬ 
tric  shadow  on  the  terrestrial  surface  with' a  plus  sign  in  the  shadow 
and  minus  in  the  illuminated  area. 

On  the  axis  of  ordinates  the  reduced  height  y  is  plotted,  and 
the  heights  in  meters  for  those  same  wavelengths  A  .  The  reflection 
formula  is  applicable  above  the  curve. 

In  drawings  36-39  the  function  V-|i(£  ,yU  ,  q)  is  given  for  q  =  o^ 
and  q=0.  On  all  the  drawings  the  variable  t  is  plotted  on  the  axis 
of  abscissas,  and  a  separate  curve  corresponds  to  each  value  of 
yt<  .  The  heavier  line  on  the  drawing  denotes  the  curves  which  corres¬ 
pond  to  f*  -  oo .  On  the  sector  0  for  q  =  o o  the  curves  are  inter¬ 

cepted  in  those  places  where,  as  investigations  show,  the  function 
Vii  is  no  longer  applicable  for  the  calculation  of  V  and  Vi  using 
formulas  (3- 15)-(3- 16)  and  (3 . l8)-(3. 19) . 

Furthermore,  in  drawings  37a  and  39a  the  actual  and  imaginary 
parts  of  function  V^^+V^)  are  given  for  q  =  oO  and  q=0  correspond¬ 
ingly.  Since  the  calculation  of  v(°)  using  formula  (3-12)  is  reduced 
to  the  use  of  the  known  tables,  then  drawings  37a  and  39a  make  it 
possible  easily  to  obtain  the  function  Vn  for  q=co  and  q=0  with 
any  values  of  the  parameter  . 

<)  V{x,  0.  v.,  oo) 

For  the  functions - -  and  V  0>^2>  makinS 

it  possible  to  find  V  for  1  using  formulas  (4.19)  and  (4.20) 

and  depending,  just  as  Vi(Z,y,q)  on  two  variables  X  and  tjg  ,  the 
3ame  graphs  are  given  as  for  Vi:  dependence  of  the  logarithms  of 
the  moduli  and  arguments  of  functions  on  X  with  fixed  values  of 
(drawings  40,  41,  46,  47);  their  dependence  on  with  fixed 
values  of  X  (drawings  42,  43,  48,  49;  and  the  level  lines  (drawings 
44,  45,  50,  51).  On  the  latter  on  the  axis  of  abscissas  the  reduced 
distance  %  is  plotted,  and  on  the  axis  of  ordinates  the  reduced 
height  . 
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On  the  following  graphs  we  give  the  functions  from  one  variable 
(with  the  same  fixed  values  of  q=©o  and  q=0).  In  drawing  52  the 

<j*  v  ix,  0,  d,  «>) 

logarithms  of  the  moduli  and  phases  of  functions  - djTdy - 

and  V ( JkT, 0,0,0)  are  given  simultaneously,  and  in  drawing  53  the  func¬ 
tions  g(z)  and  /(z),  making  it  possible  to  determine  V  for  1 

using  formulas  (4.21)  and  (4.22)  and  Vi  for<  y  4<  1  using  formulas 
(4.23)  and  (4.24). 

The  graphs  in  drawings  54  and  55  are  intended  for  determination 
of  V  and  Vi  using  monomial  formulas,  i.e.,  using  the  first  terms 
of  series  (2.6),  (2.9),  (2.7),  (2.10)  for-  large  values  of  'Stand  Z  . 

Drawing  54  gives  the  factors,  depending  on  the  reduced  distance 
*  and  Zj 
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and  in  drawing  55  -  height  factors  gi(<j)  and  /j(jj). 

Functions  V  and  Vi  are  calculated  using  the  formulas 


V  <j w)  -  a  lx,  ec)g,  (y,)g,  (y  .)  ) 
V(*»yi.y»0)— «(*».«)/iCKj)/ity*)  »’ 


v,  (z, y,  00)  -  a,  (z,  w)g,  (y) 

v  I  (z,  y,  0)  -  a,  (*.  0)/,  ( y) 


Drawing  56  gives  the  factor  Mi(y,q),  connecting  V  and  Vi  in 
the  area  of  deep  shadow,  and  also  •*  h  for  connection  of  these 
functions  between  one  another  in  the  area  of  semishadow  and  with 
function  Vn,  and  in  drawing  57  the  phase  ,  necessary  for  these 

formulas  of  connection  (3*3),  (3-7),  (3- 15)— (3* 16),  (3* 18)— (3. 19) . 
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Since  this  phase  increases  exceedingly  rapidly,  then  we  will  give 
only  its  supplement  to  the  whole  number  of  periods. 

In  order  to  facilitate  the  application  of  our  results  to  prac¬ 
tical  calculations,  we  present  here  two  schemes  in  which  it  is  indi¬ 
cated  which  drawings  and  formulas  should  be  used  with  the  different 
values  of  arguments  X  ,  ^yy ,  'Vj  or  2,  y. 


Calculation  of  V(  X  ,  ,'j£,  q) 
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Key:  (1)  Area  of  values  X  ( 2  or  Z ) ;  (2)  Area  of  values  yt  ;  (3)  Area 
of  values  Jf*  ;  (4)  Tables;  (5)  Drawings;  (6)  Formulas;  (7)  Notes; 
(8)  for;  (9)  (3-7)  applicable  in  deep  shadow  and  when  y  >  4.8; 

(10)  (3«3)  applicable  when  y  >15(10),  Z>  -2;  (11)  Initially  the 

second  pair  of  formulas  is  used,  then  the  first. 


Calculation  of  Vi  (  7  ,y  ,  q) 
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7.  Examples  of  the  Calculation  of  Field  Intensity 


We  will  give  several  examples  of  calculations  of  field  intensity 
based  on  the  tables  and  graphs  of  our  work. 

The  horizontal  component  E^>  of  field  intensity,  created  by 
a  horizontal  electric  dipole  in  the  azimuth  plane  <P  =  ~  ,  perpendic- 
ular  to  the  moment  of  the  dipole  itself,  has  the  form  [compare 
formula  (1.3)] 

£f  -  3.10* 1  !  JHtilM,  (7.1) 

where  V(  7-  ,  fa  ,  Jjt  )  is  the  attenuation  factor,  and  s  -  distance 

/^sjnir 

in  kilometers.  (The  factor  V  0  in  practice  can  be  disregarded 
within  the  limits  of  1,000  km.) 

Calculation  formula  (7.1)  gives  an  amplitude  value  of  horizontal 
component  of  the  electric  field  with  that  moment  of  the  dipole  when 
in  free  space  it  radiates  a  power  of  one  kilowatt.  If  the  power 
of  radiation  (as  this  is  often  used)  in  free  space  is  equal  to 
2  kW  and  it  is  desired  to  know  the  effective  (mean-square)  value 
of  the  field,  then  the  intensity  obtained  by  formula  (7.1)  has  to 
be  cut  in  half. 

Drawings  58-60  and  61,  62  give  the  graphs  of  the  electric  field 
according  to  formula  (7.1)  for  different  heights  of  the  correspond¬ 
ing  points  when  X  =10  cm  and  X  =7  m.  In  both  cases  the  effective 
radius  is  taken  equal  to  £L*  =  10,000  km.  On  the  left  scale  JUM/ m 
is  given,  and  on  the  right  -  the  same  values  in  decibels.  For  con¬ 
struction  of  these  curves  it  is  necessary  to  use  graphs  of  all  types: 
basic,  containing  the  attenuation  factor  V  (drawing  11  or  13),  graphs 
of  V  for  small  %  ,  and  y*  (drawing  19), 
the  functions 

dV{x. O.-Vi. oo)  V[x, 0,0, 00) 

d*3  fry  dy, 

(drawings  40-44  and  52)  for  the  case  when  one  or  both  y  are  small,  and 
Xis  not  very  small,  finally  the  functions  Vi(*,y,  ,o*)  and  Vn  (£,/«.  ,<*>), 
where  1  =  ,  and  £  =  X-Vy7 *  (drawings  24-28  and  36)  when 

44. 


one  or  both  y  are  great.  When  function  Vi  is  used  it  is  necessary 
to  observe  care  in  the  selection  of  the  transient  factor  from  Vi 
to  V,  keeping  in  mind  that  Mi  is  applicable  only  in  the  area  of  deep 
shadow,  and  ^  ^  is  applicable  in  the  areas  of  semishadow  and 

light,  however,  with  not  too  large  negative  2  s  * 

There  where  function  Vi  is  not  applicable  the  reflection  formula 
already  functions. 

Using  the  corresponding  graphs  it  is  pbssible  in  a  completely 
analogous  manner  and  in  the  same  areas  to  obtain  curves  for  the 
phase  of  the  field  components. 

Drawing  63  gives  comparative  curves  of  field  intensity  for 
values  of  effective  radius  cl- &=6 , 370  km  and  ^s4<t=25,500  km,  with 
wavelengths  of  ^=3  cm,  10  m  and  100  m. 

The  examples  given  show  that  on  the  basis  of  the  numerical  and 
graphic  material  collected  in  this  work  it  is  possible  for  any 
specific  case  to  calculate  the  intensity  and  phase  of  the  field, 
given  by  the  theory  of  diffraction  propagation  of  radio  waves  in 
a  homogeneous  atmosphere  above  a  spherical  Earth,  when  the  latter 
can  be  considered  an  ideal  reflector. 
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